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Abstract 

An expression for the electrical conductivity at the core of a magnetar is 
derived using Boltzmann kinetic equation with the relaxation time approxi- 
mation. The rates for the relevant scattering processes, e.g., electron-electron 
and electron-proton are evaluated in presence of strong quantizing magnetic 
fields using tree level diagrams. It is found that in presence of a strong 
quantizing magnetic field, electrical conductivity behaves like a second rank 
tensor. However, if the zeroth Landau levels are only occupied by the charged 
particles, it again behaves like a scaler of a one dimensional system. 
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1. INTRODUCTION 



With the recent observational discovery of magnetars, which are assumed to be strongly 
magnetized young neutron stars and also the possible sources of soft gamma-ray repeaters 
(SGR) and anomalous X-ray pulsars (AXP) the study on the effect of strong magnetic 

fields on dense stellar matter has got a new dimension. In the recent years a lot of work 
have been done on the effect of strong magnetic field on the equation of state of dense 
stellar matter [^,0 . Some studies have also been done on the effect of strong magnetic field 
on elementary processes (e.g., weak and electromagnetic processes) occurring at the core 
region as well as at the envelope of strongly magnetized neutron stars During the past 
few years, in a number of publications we have reported our thorough investigations on the 
effect of strong magnetic field on the stability of dense quark matter and quark-hadron phase 
transition (both first order and second order) at the core of a strongly magnetized neutron 
star P-[TT|. We have also studied the effect of strong magnetic field on dense neutron matter 
with properly modifying the a — uo — p-meson type mean-field theory. In those work we 
have developed the relativistic version of mean field theory in presence of strong quantizing 
magnetic field using Hartree and Hartree-Fock model with a — uo — p mesons exchange 



interaction 12,13 



In the recent years the transport coefficients in particular, the electrical conductivity of 
electrons at the crustal region of neutron stars have also been calculated in presence of strong 
magnetic fields P] . Unfortunately, in none of these studies the dense core region of the star 
is taken into account. Now from the observational data of SGR and AXP, the strength of 
surface magnetic field of magnetars are predicted to be > lO^^G. Then it is very easy to 
show by scalar Virial theorem that the magnetic field strength at the core region may go up 
to lO^^G. This is of course strong enough to affect most of the physical processes occurring 
at the core region. As a consequence, one of the most important physical phenomena, 
the transport properties of dense core matter, in particular, the electrical conductivity of 
electrons and transport properties of emitted neutrinos will be affected significantly by the 
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strong magnetic field. To investigate the transport properties of dense matter at tlie core of 
a neutron star, we use the standard form of Boltzmann kinetic equation with relaxation time 
approximation. Since protons are too heavy compared to electrons, we have assumed that the 
electric current is solely due to the motion of electrons. Now the relaxation time is directly 
related to the rates of the relevant (weak and electromagnetic) processes, these processes 
are strongly affected by the presence of quantizing magnetic field. As a consequence, both 
qualitative and the quantitative nature of transport coefficients, e.g., electrical conductivity, 
viscosity and heat conductivity should change significantly in presence of strong magnetic 
field. In a future publication we shall report the effect of strong magnetic field on the 



transport properties of neutrinos at the core region [14 



In the present article, we shall investigate the effect of strong magnetic field on the 
electrical conductivity at the core region of strongly magnetized neutron stars or magnetars. 
The paper is organized in the following manner: In section 2, we obtain the rates of the 
relevant processes in presence of strong magnetic field. In section 3, we shall derive an 
expression for the electrical conductivity from Boltzmann kinetic equation using relaxation 
time approximation. In the last section, we have given the conclusion and discussed the 
future perspective of the work. 



2. RATES 

To obtain electrical conductivity at the core of a magnetar, we have considered the 
following two electromagnetic processes relevant for electron transport: electron-electron 
scattering 

e(pi) + e{p2) e{ki) + e{k2) (1) 
and electron-proton scattering 

e{pi) + p{p2) e{h) + p{k2) (2) 
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where pi and ki are the initial and final four momenta. Since the rate of the electron-neutrino 
weak scattering, given by 



(3) 



is several orders of magnitude less than the electromagnetic processes, it has almost no 
significance in the electrical conductivity calculation. Of course, this is one of the most 
important neutrino transport processes. For the same reason, we have neglected e — e weak 
scattering. 

To obtain electrical conductivity at the core of a magnetar, we first consider e-e scattering. 
For this process, both the direct as well as exchange diagrams are considered. From the 
definition, the transition matrix element is given by 



(4) 



(5) 



where j^^ — eipf^'^ipi is the electric current produced by the scattered electron and 

is the electromagnetic field induced by the current j^{x') of the other electron and q is the 
transferred four momentum. 

Now in the case, when only the zeroth Landau levels are occupied by the electrons the 
positive energy Dirac spinor is given by 

1 
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where Ly and are the lengths along y and z directions, Eq = po = {pi + m^Y^'^ is 
the modified form of single particle energy, py and Pz are y and z components of electron 



momentum, B is the strength of magnetic field, assumed to be along positive z-direction 
(we have chosen the gauge — (0, 0, xS, 0)). The counter part of this spinor has no 
contribution when only the zeroth Landau level is occupied. Whereas, for low temperature 
case (the chemical potential ji » T, where T is the temperature of the system. In the 
numerical calculation, we have taken T — 20MeV and found that the results are almost 
temperature independent for 5MeV < T < 30MeV) the negative energy spinor solutions 
may be neglected. 

Combining eqns.(4)-(6) and after integrating over f, y', z' (which result three 5- 
functions) , we have the transition matrix element for the direct e — e scattering diagram 

Tfi = i J S{q° -Pi + kl)6{qy - piy + kiy)6{q^ - pu - hz) 

exp[iqx] expl—iq^x'] 



g2 LlLl [lQp\plklkl{p\ + m){kl + m){p'i + m){k^ + m)Yl 
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y z 

,0 uO 



exp[-i{(p2 - k2)t - {p2y - k2y)y - {p2z - k2z)z}][J„{x%^,ki[Mx)]p2,k2dx' (7) 



where 



[M^)]p,k = u{k)-f''u{p)Io;ky,kAx)Io;p^,p^{x), (8) 



u and u are respectively the positive energy spinor and the corresponding adjoint and 



Io;k„k. = — exp 



eB f kyV 



(9) 



The integrals over x and x' can very easily be evaluated by changing the variables to r and 
R, given hy r — {x — x') and R— {x + x')/2. The integration over R the gives 



TT \V2 



^ {Ply+P2yf 



2eB 
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Similarly, the integration over r is given by 

TT \l/2 



2eB 



cxp[-2eBXf,^.]Erfc[-{2eBy/'X^,,r)] (H) 



where X(^dir) = ik2y-piy-K(^dir)) /^eB , Kf^i^) = q^+ql-qo 8iridErfc{x) is the complementary 
error function. After evaluating the integrals over d'^q with the help of 5-functions we have 
the transition matrix element for direct e — e scattering 



[16p?p^A;?A;^(p? + m){p^ + m){kl + m){kl + m)]V2 
[Jn{p\,ki)\[J^,{p2,k2)]- 



'ALlLleB 



exp 



1 



2eB 

Erfc{-{2eBf/^X^^r)) 



HPly+P2y - hy - k2y)S{pu+P2z - ku - k2z) 

Which may be written as 

T)f ) = n(<^-)(27r)3<5(p; +pO - fc? - A;°)<5(pi, - k^y - k2y)5{pu+P2. - k,, - k 



(12) 



2z) 

(13) 

The transition matrix element for e — e exchange interaction can very easily be obtained 
from the direct one just by exchanging ki and k2 {ki <-> k2)- Then we have 



[16p?p^A;?A;^(p? + m)(p^ + m)(A;? + m){k^ + m)]V2 

^2)] [J2ip2, ^l)] ^j^2^2g^ 



exp 



1 



2eS 



(Ply+P22/)^ 



exp(-2eBXj,))^r/c[-(2ei?)^/2Xe.)]^^(27r) 



(ex) 



+P2 - ^1 - ^2)<^(Pl2/ +P22/ - - k2y)5{piz + P2z " A;!^ - A;22) 



(14) 



where X(ex) = (^ly — Piy — -f^(ex))/2eS and i^(ex) = ^{dir)- Which may also be written in 
the form 



Tjf = n(-)(27r)^^5(p? - fc? - k'MPiy+P2y - hy - k2y) 

S{Plz+P2z - klz - k2z) 



(15) 



Now from the elementary definition of rate per unit volume for the scattering process, we 
have the differential rate for e — e scattering 



dW 



^idrr) ^ n(-)| {2T:f5{pl+pl -kl- kl)5{j)^y+P2y - k,y - k2y)5(p,z + P2z - hz - ^^2.) 

dkiydkizdk2ydk2zdp2ydp2z 



(16) 



Hence to obtain the differential rate dW of the processes we have to calculate | n^'''^) P, 
I n*^^^) p and the cross term. Further, to compute these quantities, we have to evaluate the 
traces of the products of 7-matrices. Using the standard techniques, we have computed the 
traces of the products of 7-matrices. Then we have 

dW(dir) - pOpOj,Oj,o ^M-^iPiy + P2yf] exp(-4eBX(2^,,)) 

{Erfc[-{2eBy/^X^a,r)]V[{pi-P2)ih.k2) + (pi.A:2)(A:i.p2) - m\p2.k2) - m\p^, k^) 

+2m^]^;^5;^^<^(p? +P2 -kl- kl)5{p,y + p2y - kiy - k2y) 

Kp-i-z + Piz - kiz - k2z)dkiydkudk2ydk2zdp2ydp2z (17) 

Similarly, the exchange term is obtained by ki <-> k2 and is given by 

dW^e.) = pOpOfcOfcO exp[^(pi, +P2yf] exp(-4e^Xf,,))£;r/c[-(2e^) ^^^(e.)]^ 

[{pi-P2){ki.k2) + {Pi-ki){p2-k2) - m^{p2.ki) - w?{pi.k2) + 2m*^i 



-6{pI +pl-kl- kl)5{piy + p2y - kly - k2y)S{pu + P2z - ku - k2z) 



(27r)3 

dkiydkizdk2ydk2zdp2ydp2z (18) 
Finally the cross term whose evaluation is a bit lengthy, but straight-forward, is given 

by 



dWi 



r 1 



(cross) ; ; ^"^P 
PlP2kik2 



^ipiy+P2yT\ exp[-2eS(X2,.^) + Xj,))] 

Erfc[-{2eBy/'X(^a^r)]Erfc[-{2eBy/'X(^,^)] 

[-2{ki.k2){pi-P2) + m'^{p2-k2) + m^{p2.ki) + rn^{p2-Pi) + 

m^(pi.fci) -|- m^(pi./c2) + w?{ki.k2) — 2m^] 
1 1 

5{p\+pl -k^- k^)S{piy + p2y - kly - k2y)Sipiz + P2z - ku 



dkiydkizdk2ydk2zdp2ydp2z 



The simplified form of these rates can very easily be obtained by integrating over the 
momentum components kiz, /C22, P2y and p2z with the help of ^-functions. Then we have 
after some simple algebraic manipulation 



Erfc 



eB 



exp 



eB 



1 



(2eB)i/2' 

(P?^2 - Pizk2z){klpl - kuP2z) - m'^{plk^ - P2zk2z) - m^iplk^ - puku) + 2m^] 



(20) 



where -f^^^j^) = {piz — PfY — ipi — A*)^, A* ^-^id Pf are the electron chemical potential and 
Fermi momentum respectively, 



[2m2 + 2pV - 2puPf]^/2 
is the root of the transcendental equation + — k^ — k2 — and 



(iJ)2 



where 6 = piz +Pf 

Similarly, for the exchange interaction we have 

A r I 



dW, 



{Erfc 



exp 



-^{kly + k2yy 



exp 



kly — k2y — K(^ex) 

(2e5)V2 



n ^ 2 



-m'^{p2-ki) — m'^{pi.k2) + 2m^] 



[{Pi-P2){ki.k2) + (pi.A;i)(p2-A;2) 
1 1 



dk\ydk2y 



(If) 

k2z=kF,klz=k\,' ,k2z=klz+piz-k2z I 



\f'ik2z)\ 



where K'^^x) ~ (P'^-z ~ Pf)"^ — iPi — 1^)^ ■ Finally for the mixed term, we have 



dW, 



^"'"''^ pW2kl """"P 



-^{kly + k2yf 



exp 



2eB 



{k2y - kly - K(^dir)y 



exp 



Erfc 



2eB 



{kly — k2y — i^(ea;)) 



Erfc 



k2y kly 



(dir) 



kly k2y K(^f,x) 



(21) 



(22) 



(23) 



(2eS)V2 

-2{pi.p2){ki.k2) +m'^{p2.k2) +m^{p2.ki) + 



(2e5)i/2 

m^{pi.p2) + rn^ipi.ki) + m^{pi.k2) + m^{ki.k2) — 2rn^ 
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dhydk2y I ^^^^^^ ,fc,, = fe(«) ,k2.=k,, +Pl.-k2. l/'rfcljf _ («) ^^^^ 

In the case oi e — p scattering, given by eqn. (2) the direct diagram only contributes and the 
form of trace term will be slightly different from that of e — e direct term, given by 

32(p? + m)(A;? + m)(j9° + M*)(fc° + M*)[{pi.p2)ih.k2) + (25) 
{pi.k2){ki.p2) - m\p2, k2) - M*\pi.ki) + 2m' M*^] (26) 

where M* is the effective proton mass in the Hartree type mean field model in presence 
of strong magnetic field |]12[. In this case we have to use p^ = {pl^^ + M^)-^/^ and /c^ = 



{k\^ + M^)^/^ and make necessary changes in eqn. (20) of dW(^dir) for e — e scattering to 
obtain dW^^i})- the next section we shall use these differential rates (eqns.(20), (23), (24) 
and rfW^(rfj>)) to obtain the electrical conductivity. 



3. ELECTRICAL CONDUCTIVITY 

To compute the electrical conductivity of the magnetized stellar matter, we consider the 
simplest form of Boltzmann kinetic equation for electrons which are assumed to be slightly 
out of local thermodynamic equilibrium. The kinetic equation is then given by IjlStp^ 



§ + ^V./ + p.|^ = C (27) 

ot po dpi 



Where / is the non-equilibrium distribution function and C is the collision term, given by 
^-J (2vr)3 (27r)3 (27r)3 ^^"^^ ^^^^ + k,) 

KViy + P2y - kly - k2y)S{pu + P2z - kiz - k2z) 

I T^i I' ^ (28) 

where T = [fMl - /i)(l - /a) - /i/2(l - /3)(1 - fi)], and (1 - /,)'s are the Pauh blocking 
factors. Replacing collision term by the relaxation time approximation, we have the relevant 
portion of the kinetic equation needed to obtain the electrical conductivity. 
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where /o is the local equilibrium distribution function for electrons (Fermi distribution), 
given by 

" exp(/3(x)((p2 + + m2)V2 _ ^(^))) + i 

Here is the longitudinal momentum and p±_ = (2ueBy^^'^^ is the corresponding trans- 
verse part for the electrons, T{x)—l/j3{x) and ^{x) are the local temperature and chemical 
potential. When only the zeroth Landau levels {v — 0) are occupied by the electrons, the 
transverse momentum becomes exactly zero. 

We now consider the first order approximation, i.e., the system is very close to its lo- 
cal equilibrium configuration, then we can write f{p,x) = fo{p) + 6f{p,x) where Sf is 
a measure of its deviation from local statistical equilibrium configuration. Then we have 
^ = -(<^/){/3,o/4,o(l - /2,o) - /2,o(l - /3,o)(l " /4,o)}, whcrc /i,o's are the equilibrium dis- 
tribution function for the ith component (Fermi distribution). The we can write 

f =-/.(! (31) 
Expressing the force term in the form, pi — eEi we have after using eqn.(3) 

^ = /o(l - fo)PqEi^ (32) 
r e 

Then just by inspection, it is very easy to realize that the relaxation time is obtained by 
evaluating the integrals over dkiy and dk2y in the expressions for the rates and the relation 
with the total rate of the processes is given by 

r{Pu) = ^ (33) 

where W = W^g + W^gJ + + W^g therefore, 

11111 



(ee) (ee) (ee) (ep) 

(dir) {ex) {mix) {dir) 
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the well know Matthiessen's rule |17]. Now the expression for electric current is given by 



J. = e//(a:,p)^A^ (35) 

Because of randomness in the system, the local current vanishes at equilibrium. Then the 
nonzero part of electric current is obtained by using the first order deviation 6f{p, x) and is 
given by 

J. = ^/fr/o(l-/o)|rf^j>i?, (36) 

In presence of a strong magnetic field of strength B > the quantum limit of magnetic 
field, where B^/^ = 4.4 x lO^^G, we have d^p = 2neBdpz. Then writing Ohm's law in the 
form ji = (TijEj, we have 

e'^(3eB 



27r2 



Jrfo{l-fo)^dpz (37) 



Which is second rank tensor with the components a^z, (^zLi <^Lz and cr_|_^. Since we are 
interested for v = Q case only, i.e., when the lowest Landau levels are populated, we have 
= cr±^ = = and 



27r' 



-eB J r{pz)fo{l-fo)^dpz (38) 



which is the only non-zero component. The electron transport in this special case becomes 
essentially one dimensional- along the direction of magnetic field. 

In fig.(l) we have plotted the variation of with electron density for (a)B = 10^ x Bj^^\ 
(h)B = 5 X 10^ X {c)B = 10^ X 5^^), {d)B = 5 x 10^ x S^*^) and (e) 5 = 10^ x B^/\ where 
j^i<^) ~ 4.4 X the quantum limit for electrons. In fig. (2) we have shown the variation 

of the nonzero component of electrical conductivity azz with magnetic field strength for 
{a)ne = n^, (b)ne = lOn^, {c)ne = 50n^, and {d)ne = lOOrig, where = 10~^fm~^, the 
typical electron density in a non-magnetic neutron star. In both these cases we have not 
taken f] equilibrium condition into account. The density of electron is in some sense arbitrary. 
Therefore, we now consider an interacting n — p — e system in /3-equilibrium, which is the 
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real physical picture at the core region of a neutron star. Then we have Up — Ug, the charge 
neutrality condition, /i„ = fip+fig, the /9-equilibrium condition (we assume that the neutrinos 
are non-degenerate, leave the system as soon as they are produced) and baryon number 
density ub — rip + Un remains invariant in electromagnetic and weak interaction processes. 
Solving these constraints self-consistently at the core of a magnetar, we obtain the chemical 
potential of the constituents and hence the density of electrons for various baryon densities 
and magnetic field strengths. We have assumed a fixed temperature T — 20MeV. In fig. (3) 
we have shown the variation of a^z with the strength of magnetic field for a fixed baryon 
number density (in the numerical computation we have taken ub = 4no, no = 0.17fm~^) 
assuming that the system is in /^-equilibrium. In fig. (4) we have plotted a^z against the 
electron density for a fixed magnetic field strength {B — 10^ B^^^) in /^-equilibrium condition. 



4. CONCLUSION 

We have studied the electrical conductivity at the core of a magnetar using Boltzmann 
kinetic equation with the relaxation time approximation. The latter is obtained from the 
rates of the electromagnetic processes, which essentially control the electron transport in the 
medium. We have further noticed that the electrical conductivity behaves like a second rank 
tensor above the critical value (quantum limit) of magnetic field strength. However, only 
azz components is non-zero when the lowest Landau levels are occupied by the charged par- 
ticles. The system effectively becomes one-dimensional in presence of quantizing magnetic 
field- electrons can move along the field direction. The electric current vanishes in the plane 
transverse to the direction of magnetic field. Since it is impossible to compute the rates an- 
alytically (even it is very difficult to obtain numerically) with the non-zero values of Landau 
quantum numbers for all the four charged particles, we have assumed that only the zeroth 
Landau levels are occupied. This makes our like much simpler. The assumption of only 
zero Landau quantum number is also justified by the presence of intense magnetic field at 
the core region of magnetars. Fig.(l) shows that there exist cut off densities, beyond which. 
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the Pauli blocking factor supresses the electron scattering and as a result the conductivity 
becomes zero. This figure also shows that the cut off density increases with the increas of 
magnetic field strength. As the magnetic field increases, for a given baryon number density, 
the chemical potential of electrons decreases, which further reduces the Pauli suppression 
factor. Similarly, we have noticed from fig. (2) that there is a minimum value for magnetic 
field strength for a particular baryon number density below which the conductivity again 
becomes zero because of same reason as discussed above. On the other hand, in the case of 
/3-equilibrium condition, electrons are generated self-consistently. So there are no such cut 
off values for the density or magnetic field and the variation is smooth. 

With this simplified picture, we have obtained the rates and finally a^z for various baryon 
densities and magnetic field strengths. 
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FIGURES 
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He/ HeO 

FIG. 1. Variation of azz in sec"-*^ with electron density for (a)5 = 10^ x 
{h)B = 5 X 10=^ X Bf^, {c)B = 10'' x B^^, {d)B = 5 x 10^ x B^f^ and (e)B = 10^ x B^f^. 
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FIG. 2. Variation oia^z in sec with magnetic field strength, for (a)ne = n\ 
(b)ne = 10n°,(c)ne = 50n° and (d)ne = 100n°. 




10' 10' 10' 10' 

FIG. 3. Variation of 0-22 with magnetic field strength for fixed baryon density {ub = 4no). The 
system is assumed to be in /3-equilibrium. 



15 
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FIG. 4. Variation of a^z with electron density for constant magnetic field strength 
{B = 10^ X Bi""^). The syste m is assumed to be in /3-equilibrium. 
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